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TURBULENT TEMPERATURE FLUCTUATIONS AND TWO- 


DIMENS IONAL HEAT TRANSFER IN A 
UNIFORM SHEAR FLOW 
"by Jay Fox 

Lewis Research Center 


SUMMARY 

Correlation equations for statistically homogeneous fluctuations of ve- 
locity and temperature at two points in an infinite uniform shear flow are de- 
rived with allowance for a temperature gradient in an arbitrary direction in a 
plane normal to the flow direction* The initially excited Isotropic turbulence 
decays and becomes anisotropic with time. After Fourier transformations are 
introduced, the resulting spectral equations are solved for the case of weak 
turbulence wherein triple correlations are neglected compared with double cor- 
relations. Spectra of turbulent heat transfer and temperature fluctuation are 
calculated. For large nondimens ional velocity gradients, the thermal eddy dif- 
fusivity in the direction normal to the velocity gradient is much larger than 
that in the direction of the velocity gradient. The thermal eddy diffusivity 
in the velocity gradient direction tends to equal the momentum eddy diffusivity 
at large velocity gradients. 


INTRODUCTION 

Phenomenological theories of turbulence, which are reviewed in reference 1, 
have recently received support from statistical turbulence theory. In a uni- 
form shear flow with decaying turbulence, a tendency of the ratio of eddy dif- 
fusivities for heat and momentum to approach unity was found for conditions 
that correspond roughly to steady channel flow (ref. 2). Developments of this 
nature do not form a basis for supplanting phenomenological theories, which are 
the only practical means of organizing quantities of experimental evidence. 
Rather, statistical theories further the understanding of turbulence and may, 
in some instances, point the way for new extensions of the phenomenological 
theories when no experimental evidence is available. 

A uniform shear flow is described by a constant gradient of mean velocity 
in a direction normal to the flow direction. No boundaries are present. 
Transient turbulence, which is spatially homogeneous, is initially established, 
for instance, by a wire screen, and the turbulence Is later studied when it is 



weak enough for the triple correlations of velocity or temperature fluctuations 
to he neglected. 

Early statistical investigations of turbulent heat transfer were concerned 
with the isotropic turbulence that arises in the absence of a mean velocity 
gradient (refs. 3 and 4). For shear flows, numerical values of the velocity 
correlations were first presented in reference 5. Additional studies of heat 
transfer, pressure fluctuations, and velocity correlations were accomplished in 
references 2, 6, and 7. In the present effort, these studies are extended to 
include the effects of a temperature gradient with components not only in the 
direction of the velocity gradient (the subject of ref. 2) but also in the di- 
rection normal to both the velocity vector and the velocity gradient. A simi- 
lar arrangement of vectors occurs in a tube flow with circumferential varia- 
tions in heat transfer. In the following development, temperature gradient ef- 
fects are shown to be separable into components; consequently, the results of 
the present investigation supplement those of reference 2* 

Several features of stronger turbulence are present in weak turbulent 
shear flows, as shown in reference 5. Transfer between eddies of different 
sizes is present, as is production of turbulence by the action of the mean ve- 
locity gradient. The decay of velocity and temperature fluctuations proceeds, 
however, despite the production effects since they are not strong enough to 
offset dissipation effects. 
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SYMBOLS 

transverse velocity gradient, dU]_/dx 2 

dimensionless transverse velocity gradient, (t - ^o)®l/^ x 2 

transverse temperature gradients, hT/dx-g and St/5x 3, respectively 

constant that depends on initial conditions 

arbitrary points 

Prandtl number, v /cl 

instantaneous pressure 

distance from P to P 1 

distance vector from P to P T 

component of r 

average temperature 

instantaneous temperature 
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transfer term for temperature fluctuations obtained by integrating 
K^_ in eq. (25) over angular coordinates of wave number 

sphere 

time 

initial value of t 
average velocity component 
instantaneous velocity component 

fluctuating part of velocity component defined by eq. (4) 
space coordinate 
thermal diffusivity 

spectrum functions of Tug or tu^ defined by eq. (29) 

Fourier transform of tuI defined by eq, (15) 

Fourier transform of u_^t 1 defined by eq. (16) 
spectrum function of t^ defined by eq, ( 30) 

Fourier transform of tt T defined by eq* (17) 
equals 1 for i = equals 0 for i ^ j 

eddy diffusivity for momentum transfer defined by eq. (34a) 
eddy diffusivity for heat transfer defined by eqs* (34b) 

Fourier transform of Tp' defined by eq. (19) 

Fourier transform of pT ' defined by eq. (20) 
spherical coordinate in wave number space 
wave number 

dimensionless wave number^ V 1 / 2 ^ - to ) 1 / 2 K 

wave number vector 

component of wave number vector 

kinematic viscosity 

dummy variable 
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p 


density 

T fluctuating part of temperature defined in eq. (3) 

cp spherical coordinate in wave number space 

cp. . Fourier transform of u.ul defined "by eq. (18) 

ij 1 J 

Subscripts : 

i,j,k values equaling 1, 2, or 3 and designating coordinate directions 

(2), (3) scalar quantities that arise from effects of dT/dxg or BT/Bxg 
Superscripts : 

(') point P' 

( — ) average value 

(*) dimensionless quantity 


ANALYTICAL FORMJIATION 


The thermal energy equations at two points P and P' can he written for 
constant properties as 


BT 5 2 t 

cFfc <3x^ a dx^ dxjj- 


( 1 ) 


and 


5T’ _ 5 2 T' 

dt + dx k _ a 5x£ 5x-£ 


( 2 ) 


where u k and T are the instantaneous velocity component and temperature, 
respectively. Cartesian components of the position vector 5 are designated 
by the subscript k, which takes on the values 1, 2, and 3. A repeated sub- 
script on a term implies a summation of three terms corresponding to the three 
v alu es of the subscript. The symbols t and a represent time and thermal 
diffusivity, respectively. A division of Instantaneous quantities into steady 
and fluctuating components is accomplished by setting 

T = T + t (3) 


and 


u k = U k + u k 


(4) 


These relations are substituted into equation (l) , and the resulting equation 
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is averaged over a large number of systems that are macros copically the same 
hut have random fluctuating quantities that are spatially homogeneous in a 
statistical sense (ensemble average). The averaged equations are subtracted 
from the unaveraged ones with the result 


St TT St a St k ^k S 2 t 

St k dx k Sx k S:% Sx k a Sx k Sx^ 


( 5 ) 


where the bar indicates an averaged quantity. The average of a fluctuating 
component is necessarily zero. At point P ! , the equation corresponding to 
equation (5) can be visualized from equations (l) and (2). In a similar 
manner, the Navier-Stokes equations were treated in reference 5 to yield (at 
point P T ) 


Sut SU' 

St" + K +u k 







a2u j 

dx i 


( 6 ) 


An equation for Tui is obtained by multiplying equation (5) by u f . and equa- 


tion (6) by T, adding, and averaging the resulting equation. In the interest 
of brevity, the turbulence is assumed weak at this point in the analysis, so 
that the triple correlations that arise are neglected compared with the double 
correlations. None of the omitted triple correlations is different from that 
in reference 2. The averaged equation is 



Sru! 




SxT 


u k uj 


St 


+ TU,! 


SU' 


+ uv. 



1 Srp' 

p 


+ V 


S 2 tu t . 

J . „ J 
dx£ 8x£ dx k 


(7) 


where the independence of fluctuating quantities at one point from the position 
of the other point has been utilized in placing the quantities inside the 
spatial derivative signs. In homogeneous turbulence, (d/Sx^)^ = S/Sr k and 

(3/Bx^.) x ^ = -5/Sr^ if x 1 = x + r . In the present case of a single steady ve- 
locity component Uq and one velocity gradient dUq/dxg, a simplifying rela- 
tion exists : 


C°k 


StuT dUi StuT 
U k) - T~ r 2 ^ — *■ 

■ K Srjj dxg ^ SrT 


( 8 ) 


which reduces equation (7) to 
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Stu' au, Stu' 

J , J 




dUn 


sr + r 2 sr + UkU j as* + tu 2 5 u aG^ = 


1 ckrp* , v d 
p dr ^ ^ 61 ^ 

where 8-y = 1 for j = 1 and 0 for j / 1. 

A similar procedure applied to the equations corresponding to equa- 
tions (5) and (6) at P' and P yields 


( 9 ) 


chi^T 1 

dt 


au-, 


- 5t 


dUj_ cju^t ' 


+ u 2 t ' 8 ii a^ + u i u k S^T + d^ r 2 -SFJ 


l , ^^i^ 

pS?t + ^ a + v) 


(10) 


and, likewise, equation (5) at P and the corresponding equation at P' yield 

+ db^T r 2 + 1“ ( ^ fT + ™k) = 2a di* k T dr k (11) 

An additional equation is obtained by applying 5/Sxj to equation ( 6) and 
noting the continuity equation Suj/Sxj = 0. This produces 


i av .*4^ a2u K , 52u ^ 
p^psy-" dx j dx i dx j dx i 


( 12 ) 


Multiplying equation (12) by t, averaging, and introducing rj = xj - xj give 

(13) 


1 _ . 2 3tU 2 
p ^r . or . ~ dxg c3r_i 

j j 


Similarly, 


1 = 2 dIJ l du 2 T ' 

p dr^ Br^ dxg Sr^ 


(14) 


Fourier transforms are introduced: 


■nij (r) 


= y*j ( ~k) exp( ±K ’ r)d K 


(15) 


mini 


■ mill III III l ■■■ 



T^(/c)exp(i/< • r)dK 


(16) 


ui*r»(r) = f 
«/-00 

~(r) = f 
J-OQ 


- r « 

u ± ul(r) = I 9 ±j (k) 

° */*-co 


Tp'(r) = 


/• 


pT'(r) = 




6(lc)exp(il< • r)dlc 
exp(i/c • r)d/c 
£(lc)exp(il< • r)dK 
£ ' (k) exp( lie • r ) die 




(17) 

(18) 

(19) 

( 20 ) 


t — i 

■where | yc| = k is the wave number, which can he interpreted as the reciprocal 
of an eddy size. 

The Fourier transforms of equations (9) and (13) are 

Sr i auq 5ri St du i l 

'ST " K 1 5*2 + S* + 5 U r 2 = - p ±K ^ - < a + v ^ K T 2 < 21 ) 


and 


i * 1*1 

-- p ^ = 2 ~T-r 2 ^- 


(22) 


where dU]_/dx 2 c)T /Sx^ are constants. Equation (22) can be subtracted 

from equation (21) and £ can thus be eliminated from the equations. A simi- 
lar procedure applied to equations (10) and (14) yields an equation of the same 
form for In the present study, the case of r = 0, ST/dxg 0, and 

<$T/cbc 3 £ 0 is considered so that X\_ = Tj. The symbols a = dU^/dxg, 
b = ST/Sxg, and c = St/cSxj are introduced so that the final equations become 


cir 2 Sr 2 
w - aK i^ 


Of n 


bcp 22 - ccp 23 + |2a 2 

K. 




K i K z 


3 Sr 3 k i k< 5 

St“ " aK l "Sk£ = " ^32 " CCp 33 + 2a —3“ r 2 


K 


(l? + 1 )>' k2 ] 

- 


vk 2 t 3 


(23) 

(24) 


and 
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( 25 ) 


st " aK l 1% = " 2br 2 “ 2c-)r 3 " 2a/c2s 
where equation (25) is the Fourier transform of equation (ll). 

SOLUTION OF SPECTRAL EQUATIONS 

Isotropic turbulence and zero temperature fluctuations are assumed as 
initial conditions. Expressions for cp 22 and cp 33 that satisfy these initial 
conditions have been reported in references 5 and 7. The latter is 



where Jo and to are constants that depend on the initial conditions. 

In reference 7, the so lution for cp 2 3 = 932? although nonzero, was found 
to produce a zero value of u 2 u 3 , which was consistent with the lack of a ve- 
locity gradient in the xg^xg-plane. Equations (23) and (24) of this investiga- 
tion have been solved for the effects of cpg3 hy omission of the terms con- 
taining cp 22 and cp 33 . (The linearity of the e quati ons permits the addition 
of solutions.) Zero contributions to rug and TU3 are obtained from the di- 
rect effects of CP23; however, an indirect effect that does contribute to tu 3 
enters equation (24) in the fifth term. This contribution can be traced to the 
expression of the pressure effect {; in terms of in equation (22). The 

remaining portion of Xz that contributes to TU2 is the same as that re- 
ported in reference 2; it is not repeated herein. In the following solutions 
to equations (24) and (25), only those expressions that contribute to 'tu 3 or 
t 2 are shown. 

For a Prandtl number equal to 1 the Fourier transform of tu 3 is 
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r 3 = 


J o c + [f 


<2 + aK l( t - t o ) 
12rt 2 a k-^k 2 + k|^ 


O 2 * 4? 


exp|-2v(t - t 0 ) j K 2 + aK ± K 2 (t - t 0 ) + i a 2 /cf(t - t 0 )^j| 


x <- 


&*-,_( t - t 0 ) 


2 2 
*3*2 


<2 + |k 2 + a/c 1 (t - t Q )] 2 + /c| ( K i + ^l) 1 ^ 2 ^ 2 


r -1 k 2 . *2 + a^l(t - t 0 )1 

71 ~ 71 — i772 

L (*! + 4) (4 + 4) J 


4 K 2 


{4 + 4 ) 


tan -1 

( K l 


K 2 


+ 4) 


TJ2 


- tan 


_ 1 k 2 + a*i(t - t 0 )~| 


(4 - i) 1/2 


For a Prandtl number not equal to 1, Y "3 takes the form 


(27a) 


J 0 c{ki + [ 52 + a/c l( t - t 0 )] 2 + *|} 

12j7al^c|77c|^ 

J[ ( l/Pr ) - l]v/c 2 


expS 


aK l 


( K 1 + t- + *i) 


{ 


- 2v(t - to) [k 2 + a/c 1 /c 2 (t - to) + i a 2 /c 2 (t - t Q )^j 

»/c 2 +a/c 1 (t-t 0 ) 


exp-< 


1^1 + [*2 + a/c l( t " “to)] 2 + *§ 

(*! + ^i ) 1/2 


+ # 


s 

’ * 2 + £ 2 + 


|tan“ 


(4 + 4) 


172 


- tan' 


-1 k 2 


+ a^t - t Q ) \~l 

(*1 + K z ) lj>2 /-J 


[^ 2 + p + 1/2 ( tan (/C§ + K 2 ^ 1 / 2 


tan 


_! k 2 + aKx(t - t 0 ) N 

(4 + *l ) 1/2 




(27b) 
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These expressions for fg V3 verify the fact that the turbulent 

heat-transfer components tu 2 and tu 3 arise from temperature gradient com- 
ponents in the respective directions dT/dxg and 

The transform of the portion of T 2 that arises from St/3x 3 is, for a 
Prandtl number equal to 1, 


+ K% 


’(3) 


J o c2 { 7c 1 + [ k 2 + a /c l( t " 2 

12jt 2 a 2 k \{ k ^ + k|) 

X exp ^ - 2v(t - ■to) [* 2 + aK l K 2( t ~ t 0 ) + | a 2 K 2 (t 
a 2 Xi(t - t Q ) 2 


“ t Q ) 


5} 


+1^2 + a/c 1 (t - t Q )J 2 + 


+ ^5^1 

4 K 1 + *•! 


X 


tan 


-1 


(4 + 4 ) 


172 


. -1 k 2 

- tan x — 


+ a/c 1 (t - t 0 )T 

+ *!) 1/2 J 


(28) 


The other portion ^(2) "that is a result of ST/chcg is the same as that re- 
ported in reference 2. 


The convenience of interpreting functions of k instead of functions of 
K was pointed out in reference 8. Following that suggestion^ the integrations 
that lead to Tu 3 are accomplished in two steps, the first of which 

involves integrating over the angular coordinates of a wave number sphere: 


and 


r 3 (*) 


A( 3) ( k ) 



r 3 ( KT^CP, ©) AC 2 


sin 0 d9 d0 


5( 3 ) (/c,cp, 0)/c 2 sin 0 dcp d0 


(29) 


(30) 


A display of T 3 and A(3) shows how the contributions to 
are distributed among wave numbers or eddy sizes , since 


vu^ and 


T 


2 

(3) 


tu 3 = 


/ oo 

r 



(31) 
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COMPUTED SPECTRA 


Numerically calculated spectra of tu^ and t? . are displayed in di- 

mensionless form in figures 1 and 2 for several values of the dimensionless ve- 
locity gradient a** Since time enters all the dimensionless representations, 



0 12 3 4 


K * =, 1/2 <t-t 0 )l/2,C 

(a) Prandtl number, 1. 

Figure 1. - Dimensionless spectra of TU 3 and for uniform transverse 
velocity and temperature gradients and for various Prandtl numbers. 

the curves for various a* show the effect of velocity gradient on the spectra 
at any given time while the turbulence decays. Dashed curves correspond to 
those in reference 2 because of separability of solutions. For large Prandtl 
numbers, the spectra of TU3 in figure 1 peak at large wave numbers (small 
eddy sizes). 

Isotropic spectra (a* = 0) in figure 1 are the same for tu 3 and ru ^ 
as previously reported in reference 4. The behavior of the peaks of the spec- 
tra of TU3 and ru^ is similar to that of the respective production terms 
CCP33 and bcp 2 2 in equations (23) and (24) j 9 2 2 decreases and shifts toward 
lower wave numbers as the velocity gradient increases, whereas 933 increases 

11 


I 



markedly with little shift (see fig. 5 of ref. 5 and fig. 2 of ref. 7). 

A shift to higher wave numbers in the spectra of with increasing ve- 

- locity gradient is evident in figure 2 hoth at peak and at moderate values on 



(b) Prandtl number, 0.0L 

Figure 1. - Continued. Dimensionless spectra of ru^ and rU 2 for uni- 
form transverse velocity and temperature gradients and for various Prandtl 
numbers. 

the high wave number side; the shift results in an elongation of the spectra to- 
ward high wave numbers. This spectral change is evidently due to a transfer of 
activity from low wave numbers (large eddies) to high wave numbers (small 
eddies) by the action of the second term in equation (25), which is known as 
the transfer term. The name stems from the Fourier transform relation 
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(32) 


r 2 


Stt' 

Sri 



Ki 


5s 

S/C 2 


exp(iAT 


r)dK 


which becomes, for r = 0, 


1 




d X = 


0 


(33) 


Similar results can "be obtained from corresponding te rms in e quati ons (23) 
and (24). Thus, these terms contribute nothing to Sru^St, cHu 2 /St, and 

St^/St, but they do alter spectral distributions. 


The integration shown in equation (33) can be accomplished in two steps by 



(c) Prandtl number, 10. 


Figure 1. - Concluded. Dimensionless spectra of tu 3 and r\s^ for uni- 
form transverse velocity and temperature gradients and for various Prandtl 
numbers. 
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Figure 2. - Dimensionless spectra of l 2 an dT 2 for uniform transverse veiocity and temperature 
gradients. Prandtl number, 1. ' t2) 


first integrating over the angular coordinates of a wave number sphere and then 
integrating over the wave numbers. Of course, the second step yields a trivial 
result, hut the first result is a spectral transfer function. The integrated 
transfer term corresponding to t 2 is shown in figure 3 for a Prandtl number 
of 1. Most of the transfer of activity is out of the low wave number spectrum 
and into the high wave number spectrum, but some reverse transfer occurs at low 
wave numbers and low velocity gradients. In reference 2 this activity transfer 
was attributed to a vortex-stretching process, which might also involve vortex 
compression at low velocity gradients and thereby produce some reverse trans- 
fer. 


PRODUCTION, TEMPERATURE FLUCTUATION, 

AND CONDUCTION SPECTRA 

Production of temperature fluctuations by the third and the fourth terms 
in equation (25) is interpreted as a result of the action of the temperature 
gradient on the respective turbulent heat transfer, TU2 and TU3. Conduction 
or dissipation in the last term reduces local temperature peaks by molecular 
heat conduction away from hot spots. Production and conduction terms can be 

integrated over a wave number sphere to yield spectral distributions in the 

same manner as that used to obtain the temperature fluctuation spectra of t 2 
in figure 2. After normalization of the peak values to unity, all three spec- 
tra are shown in figure 4 for a Prandtl number of 1 and a high velocity gradient 
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a* - It -t 0 ) dU x /dx 2 


\ / 

\ / 
Vy 


"Col c Z' T+f J nC 2 ) 


tJ( 2 ) *^ 2 Tt( 2 )/U 0 b 2 ) from ref. [2] 


I I I I I I I ._l I 

.4 .8 1.2 1.6 2.0 2.4 2.8 3.2 3.6 

K * “!/ 1/2 (t -t Q ) ll2 K 

Figure 3. - Dimensionless spectra of transfer terms in spectral equations for T 2 and T 2 . Pranc 
number, 1. . (3) (2) 


^Fluctuation 


t Production ^Dissipation 

.4 .8 1.2 1.6 2.0 2.4 2.8 3.2 3.6 

k * -i/ 1/2 (t -t 0 ) 1/2 k 

Figure 4. - Comparison of production, temperature fluctuation, and conduction spectra from spectral 

equations for T 2 ^ and T 2 ^ (solid and dashed curves, respectively) at large velocity gradient. 
Prandtl number, 1. a* * (t - t Q ) dl^^ = 50. (Curves normalized to same height.) 


x 



(a* = 50) . Actually, two sets of spectra corresponding to the separate effects 
of Bt/5x 3 and ^T/oxg (from ref* 2) are displayed in figure 4. For low ve- 
locity gradients, the three spectra are close together, like those in fig- 
ure 4 of reference 2. For a large velocity gradient, production, fluctuation, 
and conduction spectra in the present figure 4 peak at successively greater 
wave numbers* 

All these effects take place as the turbulence and the turbulent tempera- 
ture fluctuations decay. Fluctuations are produced in the large eddies (low 
wave numbers), transferred to the small eddies (high wave numbers), and finally 
dissipated by molecular conduction. 



Figure 5. - Temperature-velocity correlation coefficients as a function of dimensionless velocity gradient. 
Prandtl number, 1. 


TEMPERATURE - VELOCITY CORRELATION COEFFICIENT 

A temperature -velocity correlation coefficient introduced in reference 3 
is modified herein to account for the separate effects of ST/Sxg and 5 t/Sx 3* 
Two dimensionless coefficients are utilized, 

^2 

("Lf 

and 



16 



the former being the same as that presented in reference 2. The latter coef- 
ficient has been calculated from integrals of the curves in figures 1 and 2 and 
those in figure 2 of reference 7 , all for a Erandtl number of 1. Figure 5 is 
a display of both correlations as a function of velocity gradient, starting 
with the perfect correlation value of -1 that was obtained in reference 4 for 
isotropic turbulence (a* = 0). One correlation coefficient 


tu 3 



achieves an asymptotic value of -0.9, whereas the other, by decreasing mono- 
tonically, shows a continuous loss of correlation between the temperature and 
velocity fluctuations as a* increases. 



Figure 6. - Ratio of eddy diffusivity for heat transfer to that for momentum transfer 
as a function of dimensionless velocity gradient. 


EDDY DIFFUSIVITIES 

The eddy diffusivities of momentum and heat (in the x->- and X3-directions) 
are defined hy 
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€ 


(34a) 



(34b) 


Ratios of eddy dif fusivities play a large part in phenomenological theories of 
steady turbulent flows, A unity value of %(2)/ € produces the best agreement 
between experiment and analysis for Prandtl numbers that are not too low 
(ref. 1). In the transient turbulence analysis of reference 2, a similar ten- 
dency of %(2)/ e toward unity is obtained at high values of a*, which cor- 
respond roughly to steady turbulent flows. Recent phenomenological analyses 
(refs. 9 and 10) of circumferential variations of heat transfer in round tubes 
are based on an assumption of equal eddy diffusivities in the radial and the 
circumferential directions; that is, e h(2) = %(3) in ’ tiie P r ^sent notation, 

A dimensionless eddy diffusivity v 5 A*(t - t 0 ) 3 / 2 e h(3)/ J o can be ob- 
tained by integrating the curves In figure 1. Integration of the curves in 
reference 5 for € is also necessary for the calculation of %(3)/ a ; which Is 
displayed in figure 6 along with € h(2)/ 6 from reference 2. Although the 
curves for the two ratios are not widely separated at low velocity gradients, 
which are near the isotropic case (a* = 0), large velocity gradients produce 
values of e h(3)/ € ’fcha't ara two orders of magnitude greater than values of 
e h(2)/ € ; except for low Prandtl numbers. 

The relative magnitudes of %(3) end %(2) can compared with the 
magnitudes of the turbulent velocity fluctuations (or turbulent energy compo- 
nents) in the two directions u| and u|. References 5 and 7 show that u|" 

proceeds rapidly but asymptotically toward zero at large velocity gradients, 
whereas u^ decreases slowly from the average of the energy components 
u i u i/3j "wM-ch. Increases with velocity gradient. Likewise , from physical rea- 
soning, it is clear that the thermal eddy diffusivity is greater in the direc- 
tion of greater velocity fluctuations. 

The existence of greater u| than u| has long been suspected (ref. ll) 

and. In recent times, has been verified experimentally in tube and channel flow 
(refs. 12 and 13) and In boundary layers (ref. 14). In fact, the ordering of 
all three components of turbulent energy, from largest to smallest, is the same 

^u|, u u|^ in those measurements and in the present theory (ref. 7) at large 

velocity gradients. Apparently, not all features of boundary layers and tube 
flow depend on the presence of boundaries. 
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The transient analysis was compared in reference 2 with a steady flow in a 
boundary layer or tube "by taking K* ~ 1 from turbulent energy spectral curves 
and 0.3 6 as a representative length, where 6 is the boundary- layer thick- 
ness or the tube radius. If U is an average velocity and dUp/dx^ ~ U/b, 
then a* is of the order of 0.1 Ub/v. This implies that e h(3)/ e rnuc ^ L 
larger than e h(2)/ e for Reynolds numbers of 10^ and over that are encountered 
in practice. 

The results of the present analysis, together with existing velocity- 
fluctuation measurements, provide no support for an assumption of equal thermal 
eddy diffusivities in the radial and the circumferential directions 
( € h(2) = € h( 3) ) in turbulent tube flow. Instead, the relation € h(3) > e h(2) 
is indicated. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, February 27, 1964 
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